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The kinetics of crystallization of a binary mixture of n-paramns is treated using a theory for the 
rate of growth of chains in multicomponent systems. The kinetic chain is considered to be a strip of 
crystalline material composed of molecules of both components that is growing on a substrate of uni- 
form thickness. This substrate is considered to be a close-packed surface step. Using the binary 
mixture of C24H50 and C 2 «H 5 4 as an example, the calculations are carried out by an iterative procedure 
so that the mean thickness of the strip equals that of the substrate. This procedure necessitates modi- 
fication of the rate constants that would be used without this requirement. The rate of growth of the 
strips (total flux) is calculated as well as their composition over an appropriate range of temperature 
and a complete range of liquid compositions. The pair distribution in the strips is also determined. 
The substrate is assumed to be unkinked in the calculations, whereas the substrate is probably highly 
kinked in the actual system. For this reason the calculated total flux is probably not representative of 
the actual crystal growth rate in a binary mixture of C24H50 and C26H54. It is believed that the cal- 
culated compositions are reasonably accurate. For a binary mixture of longer rc-paraffins (e.g., C50H102 
C52H106) the substrate should have many fewer kinks, and the total flux calculated by this theory will 
be more closely related to actual crystal growth rate studies. However, no data exist for such systems. 

Key Words: Crystal growth, multicomponent systems, kinetics, rc-parafnns, phase diagrams, 
pair distributions. 



1. Introduction 

In a recent paper (subsequently abbreviated LDP) 
[l] 1 we presented a theory by which the rate of growth 
and composition of a chain growing in a multicom- 
ponent system may be calculated. There are several 
kinetic problems which may be treated by the method 
developed in LDP. In particular the chain may be 
regarded as a layer (or strip) of atoms or molecules 
at a step on a growing crystal surface, the beginning 
and end of the chain, or strip, being kinks. The 
initiation and growth of such strips is one of the 
processes operative in crystal growth. In fact, in 
theories of polymer crystallization [2,3,4] it is gen- 
erally considered that formation of such strips is the 
rate controlling step in the growth of polymer crystals, 
and hence the calculation of their rate of growth 
amounts to the calculation of the rate of growth of 
these crystals. In polymer crystals the growing step 
(or lamellar edge) is not expected to be very highly 
kinked, and this approach appears valid. We treat 
this problem in the following paper. In atomic systems, 
the situation is more complicated [5]. The rate- 
controlling step is almost certainly not the initiation 
and growth of these strips, since the step is highly 



1 Figures in brackets indicate the literature references at the end of this paper. 



kinked [5], and the problem of the calculation of the 
rate of growth of a crystal becomes essentially a 
diffusion problem. In systems of intermediate com- 
plexity, such as those treated here, the relevance of 
the kinetics of formation of these strips to the problem 
of crystal growth is somewhat harder to assess. 
However, in multicomponent systems, such as we are 
treating, it seems relatively clear that the elementary 
processes we treat in calculating the rate of growth 
and composition of these chains must be important in 
determining the composition of the crystal, for it 
must in large measure be the neighbor interaction 
energies that determine the composition. We thus 
present an application of the LDP theory to the 
problem of crystal growth of binary mixtures of 
n-paraffins. The pure normal paraffins are known to 
grow, at least from solution, by a screw-dislocation 
mechanism [6], and hence a surface step is present. 
It is expected that the composition of the chains 
calculated by this method will be a more accurate 
representation of the composition of the crystal than 
the calculated kinetics will be of the rate of crystal 
growth. However, it is expected that as the length of 
paraffin molecule increases, the calculated kinetics 
become a more adequate representation of the rate 
of crystal growth, for the growing step becomes less 
and less kinked. 
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Although concerned with a specific system, our aim 
will be less to account for the behavior of this specific 
system than to illustrate the general behavior of binary 
systems insofar as this can be done by the LDP 
theory. Considerable accurate data on the phase- 
diagrams of Ai-paraffin systems are available [7, 8] 
although, to our knowledge, no data are available on 
the kinetics of crystal growth in such mixtures. The 
homogeneous nucleation rate has been studied in 
pure ra-parafnns [9], but not the crystal growth rate. 
However, many of the energy considerations relevant 
to crystal growth are well known [10]. The solid- 
phase behavior of the ^-paraffins is also well known 
[10], and, since the LDP theory is applicable only to 
solid solutions, we will concentrate on a system in 
which solid-solutions are found at all concentrations. 
In particular, we will be concerned with the solidi- 
fication to the hexagonal rotator phase in the system 
C24H50 — C26H54, which meets these requirements. 
We will ignore any solid-solid phase transitions, and 
our task will be to calculate the rate of growth of the 
strips discussed above and their composition as a 
function of temperature and liquid composition. In 
what follows we use the words "chain" and "strip" 
interchangeably. 

In order to describe the method of calculation, we 
will have to describe those aspects of the LDP theory 
necessary to make the subsequent exposition clear. 
The LDP theory begins by assuming only nearest 
neighbor interactions, with rate constants a ij and 
/3 ij , for, respectively, adding species j to a chain ending 
in species i, and removing species j from a chain 
whose terminal elements are originally i and j. These 
rate constants are assumed to be independent of the 
length of the chain. However, the rate constants for the 
initiation of the chain (i.e., the first step in the chain) 
are different, and are called a/ } and /3-j. That is, a ( { 
gives the rate of initiation of new chains from the 
melt or solution, and /3[ gives the rate of dissolution 
of this first step back to melt or solution. 

Although there is no limit to the number of com- 
ponents or order of distribution which may be handled 
by the LDP theory, in this paper we will be con- 
cerned only with binary systems and pair distributions. 
Thus, we define N J V and P'l'' as being, respectively, 
the number of chains in the system v units long with 
species j in the v\\\ position, and i and j in ( v — l)th and 
v\S\ positions. In this application, of course, i and j 
can take on only the values 1 and 2, with 1 repre- 
senting C26H54 and 2 representing C24H50. Clearly, 
we must have 



Nl'- 



! 2^ 



(i) 



bers." Kinetic equations for the growth of chains may 
be formulated [1] in terms of these and more complex 
quantities. These may be solved at steady state under 
certain fairly broad conditions. For this purpose we 
define fluxes SI and S'J as being the net rate of growth 
of chains (at steady state) from v — 1 to v units by, re- 
spectively, adding species j to any chain, and to a chain 
whose terminal (^ — l)th element is i. We obviously 
have 



S .,=2 S ,=2S(;, 



(3) 



where St is the total flux. Thus, we may consider that 
the flux Sj, is associated with the occupation number 
N{, and the flux S l J with PV. 

Under certain and quite unrestrictive conditions on 
the rate constants, it is shown by LDP that the ratio 
of flux to its associated occupation number is inde- 
pendent of chain length (for sufficiently long chains) 
and dependent only on the terminal species in the 
chain. That is, 



Si 
NL 



~P'J 



= \j 



(4) 



where the ratio is denoted XJ. This in turn is given by 
the equations 



ryi'X' 

jLi OjI. 



<?p»+\.< 

Moreover, it is also shown in LDP that 



(5) 



#i=2 



a'-iNi, 



P'j + \J 



(6) 



For sufficiently large v, this recursion relation is 
shown to become [1] 



N->-- 






(7) 



Thus for sufficiently long chains, the occupation num- 
bers become independent of chain length, and be- 
cause of eq (4) the flux also does. "Sufficiently long" 
in this context is of the order of ten, and thus if the 
average number of lattice positions between kinks is 
of this order, this equation will be valid. 

Equation (6) is a recursion relation for the N v . When 
combined with the equation for the first step in the 
chain 



Ni^ty^py, 



(2) 



where Nl is the total number of chains v units long. 
We call these quantities, NL and P'J, "occupation num- 



Ni = 



Pi+\J 



(8) 



it is entirely possible to build up the chain on a high- 
speed computer and thus get the initial ^-dependent 
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solution as well as the ^-independent limit. (In eq (8) 
OjjTV^j is the rate of initiation of new chains from the 
melt or solution.) 

Now, having the v independent limit of N j either by 
direct solution of eq (7) or iteration of eq (8), the flux 
is easily obtainable from eq (4): 



Si = XJNJ 
S r =^kJW 



(9) 
(10) 



2. The Rate Constants 



2.1. Preliminary Considerations 

It is shown by LDP that the ratio of the rate con- 
stants is given by 



aij u 



-ev)lkT\ 



(16) 



and 



It is shown by LDP that the total flux is conserved. 
Hence, we also have 



St = £ ^W j i= 2 



ft + XJ 



(11) 



In general we will not know the exact numerical value 
of aiJV j {) or p[ and hence will not be able to predict the 
magnitude of the flux. We hope, however, to be able 
to predict correctly its dependence on temperature. 
Now, the fraction of the y'th species in the chain, fy 
is given by 



fy 



S>_ 
St 



kW 



(12) 



and is thus easily obtained. The pair distributions are 
also easily obtained. It is shown by LDP that 



p u = r: n *> 
" +1 0<H-W 



(13) 



and for chains that are sufficiently long that the occu- 
pation numbers become independent of v, 



Pi 



= exp[(fij-e^)lkT] 



(16a) 



where fXj is the chemical potential of the jth species, 
€jj is the change in free energy of the growing chain 
when j is added to a chain ending in species /, and 
€°. is the change in free energy for the initial element of 
the chain. There are two things to be noted about 
eq (16) and (16a). First, they were derived from con- 
sideration of an appropriate grand-canonical ensemble, 
under the principle of detailed balance and as a result 
kinetic factors such as interfacial transport terms do 
not appear. Secondly, eqs (16) and (16a) give only 
the ratio of the forward and backward rate constants 
in terms of the free energy changes involved in the 
process. Without more detailed knowledge of the 
details of the process we do not know how to apportion 
these free energy changes to the forward and back- 
ward rate constants. Rather than try to solve the 
problem of the detailed process, as has been at- 
tempted by Frank and Tosi for chain-folded polymer 
crystallization [3], we shall take a more empirical 
approach and solve the problem for various arbitrary 
apportionments. At equilibrium the results will, of 
course, be independent of this, and our hope is that 
this will also be the case under kinetic conditions. 
As will be seen, our hope is largely fulfilled. 



PU. 



*'" 



pv+M 



(14) 



Now, clearly the fraction fy of i, j pairs in the chain 
is given by 



fiy 



Sv = \JP i J _ 1 gVkJW 
S T ~ S T " SrP'J + XJ 



(15) 



The steps in the solution are as follows: (a) Having 
the rate constants a^ and (3 ij , compute the k-> from 
eq (5). (b) Having the k j , compute the N->, by iteration 
of eq (6) beginning with eq (8) until the Nj, arrive at a 
constant value M j . (c) Compute S j by eq (9), and 
St by eq (10). (d) Compute the composition of the 
chain by eq (12) and the pair distributions by eqs 
(14) and (15). The actual methods used for these com- 
putations will be discussed in greater detail below. 



2.2. Evaluation of fjij and e/j 

The quantity (jJLj — €ij) is the change in free energy 
of the system in removing species j from the liquid 
and placing it on the end of a chain whose terminal 
element is species £, and (jij — €?) has the analogous 
interpretation. In order to do this we need both a 
model of the liquid and a model of the growing crystal. 
Our philosophy in this is to assume the simplest model 
possible for both; considering the adequacy of the 
theory and the state of our knowledge of the kinetics 
of crystallization in these systems, to assume more 
detailed models at this state would be presumptuous. 

Accordingly, we assume that the liquid is an ideal 
solution of the two species. For the system considered 
here (C24H50 and C 2 e H54) this is not a completely 
adequate description [11] but neither is it very far 
wrong. Indeed, for our purposes, it is completely 
sufficient. We may therefore write the free energy 
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G l (T,P) of a liquid mixture of Ni molecules of C 2 6H 5 4 pure crystal at the same temperature and pressure, 
and 7V 2 molecules of C24H50 as we may write 



G l {T,P) = N 1 g\+N 2 g\+N l kT In x.+N^kT In %, (17) 



€ij = g € ; + 4acr£/(/ j - /,) + 2aaU(lj - /,) 



(19) 



where gl and gi are the free energies per molecule of 
the pure liquid components at the same temperature 
and pressure, and %\ and X2 their mole fractions. We 
thus have immediately 



Hi = g\+ kT In x\ 
H<2 = g\+ kT In x-i 



(18) 



To obtain the e/j and el we adopt a model illus- 
trated schematically in figure 1. We assume a uniform 
substrate, or crystal flake, that has already grown with 
a thickness intermediate between the extended 
lengths of the molecules, l\ and Z2. (Note that we do 
not allow any roughness in the substrate. This is dis- 
cussed more fully below.) Since the solid phase has 
an hexagonal crystal structure, we assume for sim- 
plicity that each molecule in the chain is an hexagonal 
prism as shown. We now adopt a procedure similar 
to that used by Price [4] to calculate the energies. 
Each molecule may be considered to act as a rigid 
hexagonal prism with a surface free-energy cr per 
unit area. For the difference in length under considera- 
tion, this is not a bad assumption. Then, if g € . repre- 
sents the average free energy per molecule in the 




FIGURE 1. Schematic diagram of a model for crystal growth of a 
binary mixture of n- paraffins. 

The growing strip is the kinetic chain of the LDP theory [1]. The lengths of the two 
paraffin components in extended configuration are l\ and / 2 . The cross section of the two 
components in the crystal are taken to be hexagonal since the solid solution is in the hex- 
agonal rotator phase. The strip is growing on a substrate of uniform thickness, / s . 



where 

U(x) = x x^O 

= x^O 

and a is the length of the hexagon edge. Note that we 
have not included a term for the end-surface with free 
energy cr e . This merely makes the melting point that 
of the macroscopic crystal rather than a single lamellar 
sheet and corresponds to the case of the substrate 
being a surface step. Hence we obtain for fjij — eij 



M — €ij = Agj — 4fd(TU(lj — Is) 



■2ao-U(lj-li) + kT In xj 
(20) 



where Agj is the free energy of fusion (considered as 
a positive quantity) for species j at the temperature 
in question. For the first step we obtain 



W 



e° = kgj - 4ao-f/(/j - Is) - 2acrlj + kT In xj. 



(20a) 



Equation (20) indicates that if a molecule with length 
lj > l s > U is added to a chain ending in species i, 
corrections to the free energy of fusion must be made 
because of the newly exposed surface area, but if 
lj < l s < /;, no such corrections need be made. As usual 
for small supercoolings we assume 



Ag; = 



MjA7j 
TL 



where A/ij is the heat of fusion and T J m is the melting 
point appropriate to the species in question and obtain 
finally 



flj €ij — 



AhjATj 
TL 



iaaU(lj - l s ) 



- 2arrU(lj - U) + kT In xj (21) 

where ATj=T-\ n — T, with an analogous expression for 
the first element. 

2.3. The Effect of the Substrate and the Two- 
Dimensional Character of the Crystallization Problem 

If we were interested in the strictly one-dimen- 
sional problem of the rate of growth and composition 
of a strip growing on a substrate of fixed length, then 
the free energies given by eq (20) and the ratios of 
rate constants calculated from them by eq (16) would 
be completely appropriate. It would be necessary only 
to carry out the computations as outlined at the end 
of the Introduction. This would lead to correct results 
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for a one-dimensional problem, for the LDP theory, 
within the limits of the assumptions, is rigorously 
correct. 

Crystallization, however, presents a problem of a 
two-dimensional character. What is required as a 
minimum is that the strip have the same composition 
as the substrate. In this particular case, this is the 
same as requiring that the average thickness of the 
growing strip, denned as 



where f s is the fraction of species 1 in the substrate. 
Pair or higher order distributions are neglected in the 
substrate. The average free energy of a strip of length 
v in this ensemble is given by 






l=fili+fik 



(22) 



+ l22<llr ] [eiJ-M + kTlnpl v ] (24) 



be the same as the average thickness of the substrate. 

It would appear that this problem could easily be 
solved by this procedure: (a) Choose a substrate of 
some convenient thickness, (b) Calculate, by the 
methods outlined, the composition of a strip growing on 
this substrate. Its average thickness will in general 
be different from that of the substrate, (c) Use this 
strip as a substrate for another strip, (d) Continue 
this process until the average thickness of the grown 
strip is the same as that of the substrate. This pro- 
cedure clearly gives a first and reasonable approxima- 
tion to the true state of affairs in a growing crystal. 

If this is done with the free-energies given by eq (20) 
and the rate constants calculated from them by eq (16), 
no matter how the energies are apportioned between 
the forward and backward reactions, results which 
are incorrect and aphysical are obtained. If, for ex- 
ample, a phase diagram is calculated by this formalism 
at equilibrium (total flux, Sr = 0), a minimum in the 
liquidus-composition curve occurs at a different com- 
position from that in the solidus-composition curve. 
While such a situation could be conceived off equi- 
librium (St> 0) it is clearly impossible at equilibrium. 
This situation does not occur if the computations are 
carried out at fixed substrate thickness for all compo- 
sitions of the liquid, but is caused by the varying 
thickness of the substrate and the resulting two di- 
mensional character of the problem. We must, ac- 
cordingly, modify the procedure in some manner. 

One way of proceeding would be to construct a 
two-dimensional kinetic theory. This would obviously 
be a formidable problem, and rather than attack it 
with very small probability of success, we attempt a 
more tractable problem. We seek to modify the rate- 
constant ratio (eq 16) so that at least at equilibrium 
(Sr= 0) physical results are obtained. We then assume 
that the same form of this ratio will be appropriate 
under kinetic conditions (St > 0). 

The method used for this is outlined here and de- 
veloped more fully in the appendix. An ensemble of 
strips of various lengths and various average thick- 
nesses growing on substrates of various thicknesses 
is considered. The system is considered to have mole- 
cules of two lengths, U and 4, with l\ > l>. The aver- 
age thickness of the substrate is assumed to be 



where q> v , is the fraction of strips v units long with 
species j in the first position, q[{ v is the fraction of 
strips v units long with species 7* in the -nth position 
and species i in (77 — l)th position, and p$* v is the prob- 
ability that if species i is in the (rj — l)th position in a 
strip v units long, then species 7 is in the 17th position. 
Denoting the number of chains in the ensemble by 
/YY, the total free energy is given by 



G T = 2 duN v + kT^N In (N,INi 



(25) 



The fraction / of species 1 in the strip is clearly 
given by 



'2 N v 2 </■,,, 



f= 



v=\ V=l 



(26) 



5>/v,. 



For given /Xj, €J\ and €</, the free energy given by 
eq (25) is a function of q\, tX , qV^ N V9 f H and the fraction 
/of species 1 in the strip. 

By variational methods using Lagrange multi- 
pliers, a minimum in this free-energy may be sought 
with respect to these variables, and this will give 
their values at equilibrium. This analysis is carried 
out in the appendix. 

It is shown that the proper rate constant ratios to 
obtain the correct distribution at equilibrium are 
given by 







# 



a u 



xj exp [Ahj/±TjlTi,kT - laaljjkT (2?) 



- 2{(1- 2/,+/,a>)8 1 , j- 2{/ i o>&. j] 



l,=k+MU-k) 



(23) 



— = xj exp [AhjATjITikT - 8U(j- i) 

-2Ul-2f t +f.e>)8 llJ -2tf,w8tJ (28) 
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where the quantities have the following meaning: 

C = 2aor(l ] -l 2 )/kT 

2co- 1 



Kf -x 2 ) 2j r^x 2 exp (- 
% = x x explAhtATJTlkT 



Q] 



1/2 



-A& 2 Ar 2 /ni7 t -2C(l-2/ i )] 

and Si,j is the Kronecker delta. 

At equilibrium we know that f s = f. Moreover, as 
shown in the appendix, f s = co, and under these condi- 
tions physical solutions to the problem are obtained. 
Hence, without further justification we choose the 
a y//3y given by eqs (27) and (28) to be the proper 
rate constant ratios. 

2.4. The Apportionment of the Energy 

As previously mentioned, without detailed knowl- 
edge of the microscopic processes involved, it is not 
possible to make a decision of how to apportion the 
free energy between the forward and backward re- 
actions. Since we have no such knowledge we adopt a 
much more empirical approach and carry out calcula- 
tions for various arbitrary apportionments. Thus, if 
we let the argument of the exponential in eq (28) 
be A, then we take 



a ij = Xj exp (<pA) 
piJ=exp[-(l-<p)A] 



(29a) 
(29b) 



and carry out calculations for <p = 0, 1/2 and 1. It 
will be noticed that the concentration term Xj is 
placed only in the forward reaction. This seems proper 
in that the probability of adding species j to a growing 
chain ought to be proportional to the "availability" 
of that species in the liquid. 

For the initial step, if we let the argument of the 
exponent in eq (27) be A' , we take 



<*£=*; exp 04') 



(30a) 



That is, for the initial step, we put all the energy into 
the forward reaction and take the backward reaction 
to be a constant. This affects only the magnitude of 
the flux, which, as previously mentioned, we cannot in 
any case expect to predict exactly. 

2.5. The Transport Term 

Equations (16), (27), and (28) for the rate constant 
ratios, being the result of equilibrium considera- 



tions, do not involve any transport (or diffusion) terms. 
The inclusion of such terms in eqs (29) and (30) poses 
difficult problems, for in general these terms are not 
known. If these transport terms depend upon both 
the species at the end of the chain and the species 
coming on or going off (i.e., depend upon both i and 
j) then their omission can lead to serious problems. 
Thus, if a'v and /3'^ are the rate constants including 
transport terms, the most general relation to the 
oii and fi ij are 



and 



a 'ij — a vay 



fi'ij = bijpv 



(31a) 



(31b) 



where a ij and b ij represent the transport terms. At 
equilibrium, of course, we must have a [ ^\b^ — 1. 
These new rate constants determine new \j by eq (5) 
and hence new compositions and fluxes. These co- 
efficients are, of course, unknown, and indeed in such 
theories as this are usually left as undetermined 
parameters. 

In the present case however, the similarity of the 
two species in question is such that a reasonable 
approximation to the true state of affairs may be 
obtained with a much more restrictive condition than 
eq (31). In the particular case under consideration 
here we could adopt these relations 



a'v = aav 



fi'ij = bpv g 



(32a) 



(32b) 



Now, as long as St— 0, we know the ratio a\b is unity 
from eq (16), and we note that this statement is not 
restricted to the melting temperature. At the melting 
point, of course, St=0, and a = b. For finite currents 
not too far below the melting point this latter equality 
is expected to be an excellent approximation, if not 
indeed a rigorous fact. 

Thus, taking a = b, and calculating the A/* corres- 
ponding to a'v and /3'* J we have 



(30b) whence 



A '« = y ; : 



A'* _ v oSJk'J/a 
a ~ ^jgy-h k'J/a 



(33) 



4^+ y j /c 



(34) 



and 



V* = akK 



Inspection of eq (6), (12), and (15) shows that the 
transport term as represented by the quantity a has 
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no effect on the occupation numbers and the composi- 
tion, but inspection of eq (10) shows that the flux 
calculated without the transport term must be multi- 
plied by the factor a to obtain the flux when the 
transport term is included. In the calculations carried 
out here we take a to be unity, for, as previously 
stated, we cannot expect to get the true magnitude of 
flux in any case. In addition, however, the tempera- 
ture dependence of the flux we calculate will be 
slightly wrong, and should be multiplied by an Arr- 
henius type term. For the undercoolings involved in 
this work (at most several degrees) this causes a 
negligible error in the temperature dependence. 

3. Method of Calculation 

The energies and molecular lengths chosen for the 
calculation were taken from a paper by Broadhurst 
|1()|. The lateral molecular dimensions were computed 
to give the cross-sectional area per molecule given by 
Vand [12|. The data and conditions are summarized 
in table 1. Of particular importance in the table are 
the values of surface free energy, o\ which were 
chosen to be 3, 4, and 8 ergs/cm 2 . This is lower than 
the 8—12 ergs/cm 2 obtained from the data of Turnbull 
and Cormia [9|. Since we are considering the transi- 
tion from liquid to hexagonal rotator phase we expect 
(j to be lower, and indeed, comparison with the data 
of Mazee [7] shows that it most probably is. 

The calculations were carried out on an IBM 7094 
computer as follows. 

(a) For a given value of <p and a given concentration 
#i and x-i(=\— X\) in the liquid, a starting value of 
f s — X\ was chosen and rate constants calculated by 
eqs (29), (30), (27), and (28). 

(b) From these rate constants, the X were deter- 
mined by eq (5). These were solved by an iteration 
procedure outlined in LDP. Starting with 



Bi=2««//8« 



we calculate successively 



Bi 



a l >B 



v- i 



7 £'■'+#/,. 



(35) 



(36) 



Tabu* 


1. Data an 


/ conditions 


used for calculations 


( lompound 


Length, A 


\h,lk 


T m% °K 


( loncentration 
in liquid 


C 24 H 5 « 


32.46 
35.00 


6590 
7150 


323.8 
329.5 


0-1 
0-1 



a = 2.795 \ 

cr/k\ A- 2 = 2. 172, 2.896, 5.792 

<p = 0, 1/2. 1 

Temperature: 315 °K to melting point. 



until convergence is obtained. It is shown by LDP 
that this series converges to k 1 . This is easier than 
solving eq (5) for the k j directly. Moreover, for 
n>2 such direct solution is usually impossible. 
Iteration was carried out until BJ was within 1 part in 
10 7 of B* ,. Convergence took only a few hundred 
steps except near the melting point (A. J = 0, St = 0) 
where convergence was slow (10 4 steps). The series 
does not converge for T above the melting point for 
that particular composition (X < 0, St < 0). 

(c) Having the X J , the N{ were obtained by successive 
iteration of eq (6) starting with eq (8) until convergence 
to 1 part in 10 7 was obtained. N J was chosen to be 
10 6 . 

(d) The quantities Sj, St, fj, and fy were calculated 
by eqs (9), (10), (12), and (15), respectively. 

(e) Having calculated j) we of course have the com- 
position of the grown strip. We now choose fs=f 
(the fraction of species 1 in the strip) and repeat 
the procedure, and continue this until the calculated 
/ is equal to f 8 to one part in 10 6 . 

The range of parameters given in table 1 takes about 
five minutes of computation time. 

4. Results 

What we wish to know from the computations are 

(a) The rate of growth (total flux, St) as a function 
of temperature and liquid composition. 

(b) The composition of the crystal, /, crystallizing 
from a liquid of composition x\ as a function of X\ 
and temperature. 

(c) The pair distributions /n, /12, and f?i. 

To achieve these aims, computations were carried 
out for the values of the parameters listed in table 1. 
Calculations were carried out at liquid composi- 
tions %\ from to 1 by increments of 0.1 For each 
composition calculations were carried out for all the 
integer values of temperature from 315 °K up to the 
closest integer value to the temperature at which 
Sr=0. For each of the quantities of interest we are 
particularly interested in determining the effect of 
or and (p . 

4.1. The Total Flux 

Representative results for the total flux (St) are 
given in figure 2. We have illustrated its behavior 
with temperature by giving plots of it against tempera- 
ture for cr = 4 and 8 and <p = and 1. The curves for 
the other values of the parameters are very similar. 

In all cases the total flux is essentially linear with 
temperature, but perhaps more striking is the fact 
that the slope of these curves is almost independent 
of composition. This means that the total flux is de- 
pendent essentially only on the undercooling ap- 
propriate to that composition. This is not a totally 
unexpected result, for the change in crystal free 
energy with composition (aside from entropy of mix- 
ing effects) is small considering the similarity of the 
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FIGURE 2. Plots of total flux of crystallization versus temperature in binary mixtures of 

C24H50 and C26H54. 

Each of the four figures have plots for six values of liquid concentration of C 2 6H 5 4 ranging from to 1. Figures 2(a) and 2(b) are calculated 
for o- = 4 ergs/cm 2 and for <p = and <p= 1, respectively. Figures 2(c) and 2(d) are calculated for cr=8 ergs/cm 2 and for ^^O and y=l, 
respectively. All plots are quite linear. When cr — 4 ergs/cm 2 the curves are displaced uniformly to higher temperatures with increasing 
liquid concentration of C26H54; when cr — 8 ergs/cm 2 the curves are first displaced to lower and then to higher temperatures with increasing 
liquid concentration of C26H54. This behavior results from the minimum that exists in the liquidus and solidus curves for the higher value 
of cr. 



dimensions of the molecules. To the extent that these 
curves can be interpreted as rates of crystallization, 
they imply that this rate is dependent almost entirely 
on the undercooling appropriate to that composition. 
A detailed analysis of these curves shows that their 
slopes are related systematically. If the average slope 



over the illustrated temperature interval is indicated 
by m(<p,cr) then the following approximate relations 
hold: 



m{tp, cr) 
m(0,<r) 



= /(*>) 



(37) 
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and 



m(<p,cr) 
m((p,cro) 



= g(cr) 



(38) 



where f(<p) and g{cr) are functions only of cp and cr, 
respectively. This means that the variation of these 
slopes with cp for any given cr is independent of cr 
and the variation with cr for any given <p is independent 

of (p. 

As previously discussed under the section on the 
transport terms, these slopes are not correct, for they 
do not contain the effect of this term. They should be 
multiplied by the factor a which is expected to be of 
the form 



a = aoexp [— AH/kT]> 



(39) 



However, unless the preexponential factor a and the 
activation energy AH are inordinately large, this 
will have only a negligible effect on the slopes over the 
eight degree range of undercoolings considered. 

4.2. The Crystal Composition 

The composition / of a strip crystallizing from a 
liquid with composition x\ for various values of the 
parameter cr is shown in figures 3 and 4. Figure 3 is 
computed for Si — at ail compositions and thus is a 
true phase diagram. The effect of the parameter <p 
is indistinguishable on the scale of these figures. 

Figure 3 indicates that as the surface energy of the 
material is increased, a minimum in the phase diagram 
occurs, and the maximum separation between liquidus 
and solidus increases. Since cr is a measure of the 
interaction energy between the two species in the 
solid phase, this behavior is normal. Detailed cal- 
culations indicate that the liquidus and solidus curves 
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FIGURE 3. Phase diagrams of the binary mixture 0/C24H5O and 
C2«Hs4/or three values of a. 

X\ is the molar concentration of CshH.vj. This diagram represents the composition of the 
solid at equilibrium (Sr=0) -with a given liquid composition. The solidus curve falls below 
or touches the liquidus curve. As cr increases above 4.018 ergs/cm 2 a temperature minimum 
occurs in the phase diagram. 



approach the temperature axis horizontally for a cr 
value of 4.018 ergs/cm 2 . The experimental curves 
given by Mazee [7] for this system indicate that the 
liquidus and solidus curves do not have such a hori- 
zontal portion, and we conclude that cr for the hexag- 
onal rotator phase of the normal paraffins is somewhat 
less than 4 ergs/cm 2 . This is to be compared with 
the value of about 12 ergs/cm 2 deduced from the data 
of Turnbull and Cormia [9] for the nonrotator phases. 

Figure 4 looks superficially like the phase diagram, 
figure 3, but is not, for it is for kinetic conditions 
and hence the phase rule, which is derived for equi- 
librium, cannot be applied to it. This diagram gives 
the composition of the solid separating from a given 
liquid composition at constant total flux St and we 
call it a "kinetic phase-diagram," and even use the 
terms "liquidus" and "solidus," even though not 
entirely appropriately. The constant total flux was 
chosen by choosing the crystallization temperature 
for the pure x 2 case to be 321.0 °K (undercooling of 
2.8 °C). The composition of crystal separating from a 
liquid of composition x\ was then computed for this 
same flux, as well as the temperature at which this 
occurred. By the remarks made about the behavior 
of the total flux with temperature and composition, 
this gives almost a constant amount of undercooling 
throughout the diagram, but not quite, as is apparent 
from the X\ — 1 end of the diagram. 

The curves of figure 4 are superficially very similar 
to those of figure 3, and indeed it would take extremely 
careful experiments to be able to distinguish between 
the two. This indicates that even if conditions of 
exact equilibrium do not obtain when determining 
phase diagrams for this system, very little error will 
be made. 

A more detailed analysis shows, however, that 
there are differences between these kinetic phase 
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FIGURE 4. A "Kinetic phase-diagram'' of the binary mixture of 
C24H50 and C26H54 for three values of a. 

x, is the molar concentration of C»«H H , This diagram represents the composition of the 
id separating from a given liquid concentration al constant total Mux St- The total flu\ was 
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1 temperature for the case v, = I) to be 321.0 °K (under- 
of the crystal separating from a liquid of composition 
3 well as the temperature at which this occurred. These 
a constant amount of undercooling. The effect of the 
he scale of this figure. 



253 



diagrams and the equilibrium phase diagram. The 
most striking difference is shown in figure 5 which is 
an enlargement of the region of the minimum for the 
curve calculated for cp = 1 and cr = 8 for the same total 
flux as in figure 4. This figure indicates that the 
minimum in the temperature-liquidus curve occurs 
at a different composition from that in the tempera- 
ture-solidus curve. This could not, of course, occur 
in an equilibrium phase diagram, but there is nothing 
to preclude it happening in these kinetic diagrams. 
To the right of the vertical arrow the composition of 
the crystallizing solid is higher in species one than 
that of the liquid (f>x\) while to the left it is lower 
(f<x\). This is denoted by the horizontal arrows. 
At the point denoted by the vertical arrow the com- 
position of the solid is the same as that of the liquid. 
This "equal-composition" point, which we denote by 
X e , for an equilibrium phase diagram must, of course, 
occur at the minimum of the liquidus and solidus 
curves, but no such restriction occurs here. 

It is to be noted that although the minimum in the 
liquidus and solidus curves occur at different com- 
positions, the temperatures at the minima are the 
same. For this not to be true would require a double 
valuedness in the relation between / and X\. But our 
a#, /3#, and \ j and the quantities calculated directly 
from them are continuous and single valued. Hence, 
this double-valuedness cannot exist, and the tempera- 
ture at the minimum in both the solidus and liquidus 
curves must be the same. It is also clear on physical 
grounds, that any point on the liquidus curve must 
correspond to a point on the solidus curve, and vice 
versa. Thus the minimum in the two curves must be 
at the same temperature. 

It is interesting to consider the effect of super- 
cooling on the composition X e . As the undercooling 
is increased, X e decreases, and this is illustrated 
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FIGURE 5. An enlargement about the minimum of the "kinetic 

phase diagram" of figure 4 for cr = 8 ergs/cm- and <p = l. 

Here the minima of the temperature-solidus (solid line) and the temperature-liquidus 
(dashed line) curves do not occur at the same composition. This could not occur in an 
equilibrium phase diagram. To the right of the vertical arrow the composition of the crystal- 
lizing solid is higher in C 2 6H 5 4 than that of the liquid, while to the left it is lower. This is 
indicated by the horizontal arrows. The vertical arrow indicates the "equal composition" 
point. 



schematically in figure 6. In this figure we show the 
equilibrium phase diagram and one of the family of 
kinetic phase diagrams. On this figure there are shown 
three curves intersecting the kinetic phase diagram. 
These represent schematically the locus of X e (T), 
the locus of the minimum in the liquidus curve, 
X m (T), and the locus of points on the liquidus con- 
jugate (i.e., at the same temperature) to X e . This 
last locus we denote by X c € . The form of these curves 
is difficult to develop analytically. This much, how- 
ever, can be said. The slope of X e (T) and X m (T) at 
the equilibrium diagram is negative. It does, however, 
depend upon ip and the cr involved, and this sign of 
the slope is not a general rule. The slope at the point 
where these curves cross the %i = axis is positive, 
and also depends upon <p. 




FIGURE 6. Schematic presentation of the equilibrium phase diagram 
and one of the family of "kinetic phase diagrams" of a binary 
mixture of n- paraffins. 

Three curves intersect the "kinetic phase diagram." These curves represent schemati- 
cally the locus of the equal composition point, X e (T), the locus of the minimum in the 
liquidus curve, X m (T), and the locus of the points conjugate to X {T) (i.e., at the same 
temperature on the liquidus curve), X C {T). The temperature minimum of the equilibrium 
phase diagram is T m . The supercooled area beneath the equilibrium phase diagram can be 
divided into 7 regions designated by Roman numerals and two points, a and b. In these 7 
regions and at these two points the time dependence of the total flux are qualitatively 
different. 



Moreover, the temperature T at which X e =0 can 
be determined analytically. This can be shown to be 



Afr) -A/i 2 — 6acr(/i — l 2 ) 
A/i./r.-A/i./n 



(40) 



For the parameters for the curves in figures 4 and 5, 
T= 232.5 °K which is a completely unattainable under- 
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cooling. For cr= 4.05 ergs/cm 2 and all other parameters 
the same as for the above curves, T= 322.90 which is 
an undercooling of only 0.9 deg, and hence in principle 
observable. However, in this case the minimum is 
very slight. For the system C 3 oH 6 2 — C35H72, for 
which experimental data are available and a minimum 
exists, T= 191 °K, and hence also completely un- 
attainable. We thus conclude that the existence of 
T is a prediction of the LDP theory which is likely to 
remain academic. 

The kinetic phase-diagrams in figures 4 and 5 give 
the composition of solid crystallizing at a given rate 
from a given composition of liquid. In other words, 
in a real system this gives the composition of the very 
first crystal formed. As time proceeds, however, the 
composition of the liquid changes, and it is interesting 
to deduce what happens to the total flux as time 
proceeds. 

The most interesting case is that in which the 
phase diagram has a minimum, and the behavior in 
this case can be described with reference to figures 
6 and 7. Figure 6 shows the equilibrium phase diagram 
and one of the family of constant flux kinetic phase 
diagrams. Those with smaller fluxes are above the 
illustrated kinetic phase diagram, and those with 



larger below. It is convenient to divide the tempera- 
ture-composition space into seven areas, as shown, 
and two points, a, and b. 

What happens to the total flux when a liquid of a 
given composition is brought to a temperature below 
the melting point and held at the temperature is 
shown for each of the seven regions and the two 
points in figure 7. The diagrams are almost self 
explanatory, but we will describe them briefly. 

In regions I and II, as the liquid crystallizes, its 
composition changes such that the point describing 
it moves into areas where the crystallization rate is 
lower, finally reaching a point on the equilibrium 
liquidus curve at which time the flux drops to zero 
and crystallization ceases. 

In regions IV and V, the flux initially decreases, 
but because the point representing the liquid com- 
position crosses the curve, X m (T), the flux begins to 
increase again, eventually leveling off where the 
point reaches the X € {T) curve at, respectively, lower 
and higher values. If the initial liquid composition 
and temperatures are such that the initial liquidus 
point is on the X c e {T) curve (point 6), the flux first 
drops, but then increases back to its original value 




FIGURE 7. Schematic plots of total flux, S T , versus time for a binary mixture of 
n-paraffins that is represented in figure 6. 

Each plot represents schematically the time dependence of the total flux when crystallization begins in the indicated 
regions or at the indicated points of figure 6. In regions I and II crystallization does not proceed to completion. 
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when the liquid point reaches X e (T) and remains 
constant thereafter. 

In region VI the flux increases monotonically since 
all points in this region are to the left of X m (T) even- 
tually becoming constant when the liquid point reaches 
X e (T)\ In region III, the flux decreases monotonically 
until the liquidus point reaches X e (T) after which it 
becomes constant. 

If the initial composition and temperature are such 
that the liquid point is on the curve X e (T) (point a), 
the composition and hence the flux are independent of 
time. This is clearly a generalization of the behavior 
at the minimum in the equilibrium curve. 

In region VII (all points below T) the flux decreases 
monotonically in time and reaches a constant value 
only when the species 1 is completely depleted from 
the liquid. 

It should be pointed out again that although this 
behavior is a direct consequence of the LDP treat- 
ment, its experimental verification is likely to be very 
difficult, and its experimental importance relatively 
slight. In any real system the curves X e (T), X m {T) 
and X C (T) are almost certain to be practically coinci- 
dent and vertical, and region VII will be unattainable. 
Thus, regions V, VI, and VII will practically disappear, 
and the behavior in all attainable regions will be such 
that the flux will decrease monotonically — to zero if 
the liquid point is above the minimum melting tem- 
perature but below the liquidus curve, and to a constant 
value if the liquid point is below the minimum melting 
temperature. 

The last question we have to investigate is the 
occurrence of a eutectic. It is clear from what has been 
said about the continuous single-valued nature of our 
cfi, f& and X? with respect to %\ and /, that we will 
never obtain a eutectic unless we include its possibility 
in a more or less ad hoc manner. To do this it would 
be necessary to investigate the free-energy of the solid 
phase with respect to composition and determine if 
in some composition range the free energy of the solid 
is lower when expressed as a linear combination of 
two fixed solid compositions in standard ways [13]. 
For the values of the interaction parameter (a) chosen 
here, this does not occur and no eutectic is formed. 
A simple calculation indicates that a value of cr greater 
than approximately 22 ergs/cm 2 is necessary before 
eutectic formation occurs at any reasonable under- 
cooling. This is an unreasonably high value of surface 
energy for n-paraffins. 

5. The Pair Distribution 

As indicated in the introduction the pair distribu- 
tion is easily calculated by the formalism of LDP. 
In order best to present typical results, in figure 8 we 
show plotted against composition / the difference 
between the observed /12 and the random mixing 
case, where /12 is given by /(l —f). These are com- 
puted along the solidus line in figure 4 for <r = 4 and 
8 ergs/cm 2 . This curve gives all the possible informa- 




FlGURE 8. Plot of the difference between the pair distribution 
f 11 and the random mixing pair distribution (f\) 2 versus the singlet 
distribution fi for various indicated values of cr and <p. 

These curves are calculated along the solidus curve in figure 4. The deviations from 
random mixing are very small. 

tion about the pair distributions, for the following 
identities hold: 

/l2 == /21 

/<!-/) -/ 12 =/n-/ 2 

fn-r=M-(\-ff. 

First to be noticed in figure 8 is the small deviation 
of /12 from the random distribution. In the worst case 
this amounts to only about 6 percent at the f =1 /2 
point. This near equivalence is to be expected from 
the small magnitude of the interaction energy as ex- 
pressed by 6ao-(l 1 — l 2 ) in comparison to the entropy 
of mixing terms. 

The effect of a is such as to increase the deviation 
from randomness, as is to be expected, and the 
deviation is essentially proportional to cr for a given 
value of <p. The latter parameter causes a not incon- 
siderable effect, but we have no explanation for it. 

The curves show an apparent symmetry about the 
point /= 0.5. If this symmetry were indeed present, 
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it would mean that fn is symmetric about /=0.5 
since clearly the random case is. However, this 
symmetry is only apparent. Detailed investigation 
near x = and x — 1 indicate that f 12 is not symmetric 
about /=0.5. This is what would be expected con- 
sidering the form of Equation 28, but the point is too 
recondite to pursue. 

6. Conclusions 

We have shown that by relatively routine calcula- 
tions using the LDP method, it is possible to calculate 
the kinetics of growth and composition of new layers 
growing at a surface step on a growing crystal in a 
binary mixture when the solid phase forms a solid 
solution. This procedure is expected to be an adequate 
representation of one of the processes occurring in 
crystal growth when the average distance between 
kinks is of the order of ten or more lattice spacings. 
If this distance is very large then this calculation 
would give the kinetics of crystal growth, but it is 
expected that in the present application the results 
will be more accurate for the composition. If the cal- 
culations were made for a binary mixture of longer 
hydrocarbons (e.g., C50H102 — C52H106) it may be ex- 
pected that the distance between kinks would be large, 
and the flux calculated would have relevance to the 
crystal growth rates. 

What it is necessary to know for the calculation in 
this paper are the rate constants for the addition and 
subtraction of a given species in the mixture to and 
from the growing strip. It is possible by equilibrium 
considerations to evaluate the free energy change in 
the system and relate this change to ratios of the rate 
constants for addition and subtraction, but it is not 
possible in general to evaluate the rate constants 
themselves. However, in this problem, various ap- 
portionments of this free energy change to the addition 
and subtraction rate constants have only minor effects 
on the results. 

In the process of carrying out these calculations, 
specific account had to be taken of the two-dimen- 
sional nature of crystallization. This necessitated 
adapting the LDP theory, which is an exact theory 
for one-dimensional kinetic problems, to this two- 
dimensional problem. This was accomplished in an 
approximate, but in our view adequate, manner by 
appropriate modification of the rate constants. 

Although these calculations were carried out on the 
rather simple binary system C24H50 and C26H54, there 
is in principle no reason why it could not be carried 
out on more complex mixtures, provided only that the 
interaction energy in the solid is not so great as to 
cause eutectic formation. Indeed, in the following 
paper we apply the LDP method to a system where as 
many as 50 components must be considered. 

7. Appendix 



The physical model and the concepts are emphasized, 
while many mathematical steps are omitted. 

We calculate these ratios from an ensemble that 
consists of a "mother liquor" formed by two hydro- 
carbon components that surrounds a large number of 
crystals. The crystals present Nr substrates to the 
"mother liquor" upon which strips can be formed. 
The length of a strip (the number of molecules crys- 
tallized from the "mother liquor") is designated by 
v and v can take on the values 0, 1, 2, . . ., v {) . 

We assume that all substrates have a uniform thick- 
ness, l 8 . This uniform thickness is taken as the average 
thickness of the substrate, and is given by eq (23) of 
the text. The assumption of uniform thickness is nec- 
essary in order that we can apply our two component 
one-dimensional rate theory to the problem. This 
assumption is clearly an approximation, because the 
substrates will be of uneven rather than of uniform 
height. We shall see below that this assumption causes 
a modification of the rate constants. 

The strips will be deposited upon the substrates in 
such a way that the free energy of the ensemble will 
be minimized. We shall determine the occupation 
numbers that characterize the strips such that the 
total free energy is minimized. These occupation num- 
bers can be used to evaluate the ratio {a'J/ftv) by the 
principle of detailed balance. This is analogous to the 
method used in single component rate theory. In the 
present case this leads to fl| 






m 



(1A) 



In the present instance the ratio {av/fi'J) is independent 
of v, and we have anticipated this result in writing 
eq (1A). It should be noted that the occupation num- 
bers on the right-hand side of eq (1A) are not the 
steady state occupation numbers obtained elsewhere 
in this paper. The total free energy of our ensemble is 



G r = ^NvC+kT^NAn 



Nr 



(2A) 



The reference state for this free energy expression 
is the state where no strips are deposited on the sub- 
strates. The average free energy of a strip of length v 
is given by 



G,=2 9i. d^- N + kT in ^ J 



+ E t 2 it •>!>«- M+ kT In p% ,] . (3A) 

T) = 2 i= l j= 1 J 



In this appendix we derive expressions for {a'J/l3 ij ) The quantities appearing in eqs (2A) and (3A) have 
that are suitable for our model and our treatment, all been defined in the text. The quantities appearing 
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in these equations are not independent but must 
satisfy the following relations: 



The quantity T in this equation is a Lagrange multi- 
plier associated with the variation of eq (9A). 
We also define the row vector 



^N V = N T , 



2 «*.«=i 2 2<7«„=i, 



(4A) 



(5A) 



q\, , = 2 vi'. 2 </i. » = 2 «- " +1 = 2 # 1- ( 6A ) 



In addition the conditional probabilities in Equation 
(3A) satisfy the relations 



q$,r,=qi,vPtv 



2^'.r ) =l. 



(7A) 
(8A) 



The fraction of component one in the strips is a) 
and is given by 



w' () = (w x {) , wf } ) 



(HA) 



where 



[(W.-rw+ra>] (12 A) 



(Agj 2a<rlj 

4acr(/ 1 — L) 
kT 



With these quantities We can write down our solu- 
tions as 

<(W ) M_ *(5) «*i (8i U) (F)-"/ 1 ) 

<# » = V ,^„-,/T \ -/ (13A) 



oW- 1 / 1 



?£. » : 



w 



w' a (w)^{^)(M,?>i){wy^Q 



w«{wy-iQi 



(14A) 



2^2 < 

y=j i? = i 

2 vN ' 



(9A) 



A /I 



(15A) 



These solutions allow us to evaluate the occupation 
numbers and evaluate the ratio of the rate constants. 



We discuss below the relationship between f s and a>- 
the fractions of component one in the substrate and 
strips respectively. 

We have written down the total free energy of the 
ensemble. In order to determine the occupation num- 
bers that minimize this free energy, we calculate the 
variation of 67 , allowing the N v , q{^ u qi{ rt,f s , and a> to 
vary, and equate 8Gt to zero. Since these quantities 
are not all independent, but instead are related by 
eq (4A), (5A), (6A), and (9A), we introduce Lagrange 
multipliers with these equations. Then the coefficients 
of 8iV„, 8qi p ,, and 8gj? v are equated to zero. 
This leads to quite complicated expressions, but the 
algebra is relatively straightforward. 

The procedure followed here is similar to that em- 
ployed by Mullins fl4j for Markov chains, except that 
in our case the chains have finite lengths and the 
matrices involved are not symmetric. We shall not 
write down the equations involved but merely write 
down the solutions which are most easily expressed 
in matrix notation. Let us define the W matrix 



(W)ij = Wij = Xj exp 



Agj 2aaU(l j -l i ) 



kT 



kT 






PJL \ - ^+lgfc.1.«H-1 _. 



m 



N v ql 



(16A) 



We therefore have an expression for the ratio of the 
rate constants. 

Our solution is not complete, for the elements of 
the (W) matrix contain the quantities T and co, which 
as yet have not been determined. We can obtain two 
further relationships. From eqs (9A), (14A), and (15 A) 
it can be shown that for v Q very large 



/= 



-W22 



2z — w\\ —W22 



(17A) 



where z is the largest eigenvalue of the W matrix. 
This equation provides a relationship between co 
and T. 

A further relation is obtained by considering the 
variation of the ensemble total free energy. Using 
Equations (13A), (14A), and (15A) this becomes 

8G T = 4aer ( |) vn)j (rSco - a>8f 8 ). (18A) 



~~j7f ~ [(1 — f*~ r)8/+rco] \ (10A) This equation relates the variation of the total free 
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energy of the ensemble to the variation in the composi- 
tion of the substrates and strips. 

The two dimensional character of our problem must 
be kept in mind in our treatment of this equation. 
This can be illustrated by considering a truly one di- 
mensional problem as an example: Let the mother 
liquor and strips be the two component system of 
hydrocarbons under consideration, but suppose the 
substrate consists of a third component, say still 
another hydrocarbon. Then we need consider no 
variation in f s i.e., 8f s = 0. In this case in order that 
8Gr=0, we set r=0. The ratios of the rate constants 
become the usual Boltzmann factors that were derived 
in [1], These rate constants are the proper ones for 
the one dimensional problem. However, if we follow 
this procedure in the present problem (i. e., set 
8F. S = 0, r=0), the resulting rate constants are not 
the proper ones. This is most clearly seen when the 
equilibrium phase diagram has minima in the liquidus 
and solidus curves. These minima do not occur at the 
same value of the concentration, which would imply 
a region of concentrations where no stable equilib- 
rium exists. 

We must therefore treat eq (18A) more carefully. 
Suppose we have a particular liquid concentration at 
a particular temperature. There will exist a crystal 
concentration, f<{x x J), that is most stable. That is, 
if we vary the concentration of the crystal (8f=8f s ) 



3G T = 4ao- ( f ^) ff -feW* = 

T =/„(*,, 7) = a. =/,. (19A) 



so that 



This property must hold at equilibrium (i.e., where 
x u T fall on the equilibrium liquidus curve), but it 
also will hold at other concentrations and tempera- 
tures. This result does not solve our problem, however, 
except at particular values. We wish a value of T 
that holds for any values of (x u T,f s ), while eq (19A) 
yields a value of T that holds for (#, T, f 8 =fe)- When 
the kinetics of the crystallization process leads to 
values of f s ^f e , it cannot be expected that 8G T 
vanishes. 

A more critical analysis would be required to obtain 
values of V for general values of Ui, 7\ / s ). Such an 
analysis does not seem worthwhile, because we have 
already introduced the approximation that the sub- 
strate height is uniform. Instead we shall assume 



This assumption satisfies eq (19A) when f 8 = fe, and 
should be reasonably accurate at other values of 
f 8 . At thermal equilibrium it satisfies our requirements. 
We therefore have 

aV [ A& _ 2aa-U(lj-li) 

T I 

4acr(/i — li 



pa Xjexp { 



kT 



[(1-2/^+zh}- (21A) 



In this equation a> is the value that must satisfy eq 
(17A). This can be expressed as 



2g>-1 = 



g-l-f*! 



/(g-l+x0 2 + 4(l-^gexp { 2a(7( ^ fa) } 



where 



t = f A S' A g* 4ao-(/ 1 -/ 2 )(l-2/; s .) 

Z *iexpj AT kT > kT 



(22A) 



(23A) 



r=/.. 



(20A) 



Then we see that for any values of (xi, T, f s ) we can 
obtain co from eqs (22A) and (23A). Then the ratios 
(avipv) are obtained from eq (21A). Equation (21A) 
is the same as eq (28) of the text. 
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